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I.  INTRODUCTION 

Until  recently,  few  alternatives  existed  for  the  use  of 
hydraulic  or  pneumatic  actuators  in  primary  flight  control 
applications.  However,  advances  made  in  the  field  of  rare 
earth,  permanent  magnet  materials  and  in  high  power  semicon- 
ductor transistor  technology  has  made  possible  the  use  of 
electromechanical  actuators  as  a  practical  alternative. 
Utilizing  rare  earth  (specifically,  samarium  cobalt) 
magnetic  material  within  a  brushless  dc  motor  design 
provides  a  prime  mover  in  flight  control  applications 
offering  superior  performance  characteristics  over  electroh- 
ydraulic  systems. 

Elimination  of  the  brush  type  commutation  scheme  within 
the  dc  motor  provides  numerous  advantages:  1)  higher  rated 
motor  speed  along  with  a  reduction  in  weight  and  volume  for 
a  given  horsepower,  2)  the  ability  to  use  permanent  magnet 
rotors  instead  of  a  rotating  armature  winding,  combi'ned  with 
the  elimination  of  the  brush  assembly  translates  into 
design,  implementation  and  maintenance  simplifications,  3) 
siQce  there  are  no  brushes,  no  arcing  will  occur,  hence 
allowing  motor  operation  in  hazardous  environments,  and  4) 
improved  thermal  characteristics,  as  losses  (ohmic  and 
core)  ,  which  arise  primarily  within  the  stationary  portion 
of  the  machine,  are  easily  dissipated  through  the  stator 
housing. 

Use  of  the  brushless  dc  machines  does  not  come  without 
cerain  disadvantages,  chief  among  them  being  the  cost  and 
the  uncertain  availability  of  the  samarium  cobalt  material 
for  rotor  construction.  As  an  aside,  a  recent  study  was 
conducted  in  which  the  performance  characteristics  of  both  a 
ferrite   and   a   samarium    cobalt   type    dc   motor    w*=re 
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investigated  [ Ref .  1  ]•  The  research  conducted  demonstrated 
the  superiority  of  the  samarium  cobalt  design,  and  hence  its 
desirability     for      present      and      future      applications.  An 

additional  disadvantage  accrues  from  the  fact  that  motor 
commutation  must  be  accomplished  electronically,  resulting 
in  an  increase  in  the  complexity  of  the  design  of  the  motor 
controller,  with  a  commensurate  increase  in  the  cost  of  its 
implementation. 

This  study  will  concentrate  on  the  power  control  of 
brushless  dc  motors  utilizing  a  pulsed  power  approach.  Power 
pulse  control  of  dc  motors,  better  known  as  pulsewidth  modu- 
lation (PWM),  utilizes  as  an  input  to  the  motor,  voltage  or 
current  pulses  with  the  pulse  duration  being  controlled. 
Pulsewidth  modulation  offers  considerable  advantages  in  the 
control  of  dc  motors,  which  will  be  outlined  in  the 
following  chapter.  The  research  was  conducted  utilizing  a 
computer  model  of  a  3-phase,  4  pole  brushless  dc  motor, 
which   was   developed    by  Thorn  as   in    a   related   study    [fief.    2]- 

It  is  to  be  noted  that  this  research  project  represents 
only  one  part  of  an  ongoing  effort  to  accurately  simulate 
the  use  of  a  brushless  dc  motor  as  an  electromagnetic  actu- 
ator for  use  in  advanced  missile  control  systems.  While 
electromagnetic  actuators  have  previously  seen  use  in 
missile  projects  such  as  HARM  and  Condor,  these  actuators 
have  teen  too  large  and  have  had  too  slow  a  response  for 
high  torque  applications  as  are  found  in  the  AMfiAAM 
(Advanced  Medium  Range  Air  to  Air  Missile)  project  [Ref.  3]. 
The  overall  aim  of  the  entire  research  program  is  to  help 
exploit  the  technological  opportunity  which  exists  within 
rare  earth  permanent  magnet  dc  machines  in  the  role  of  elec- 
tromechanical actuators  for  use  in  guided  missile  control 
appl ications. 
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II.     PULSEWIDTH    MODULATION 

Pulsewidth  modulated  switching  amplifiers  offer  consid- 
erable advantages  in  the  control  of  dc  motors.  Before 
describing  the  techniques  used  in  simulating  the  pulsewidth 
modulated  control  of  the  computer  model  mentioned  previ- 
ously/ it  will  be  important  first  to  outline  the  basic  prin- 
ciples  and    characteristics    of    pulsewidth   modulation. 

A.       PULSEWIDTH    MODULATION    PRINCIPLES 

The  pulsewidth  modulation  scheme  utilizes  transistors  in 
the  switching  mode,  whereby  the  transistors  are  switched 
into    and   out  of    saturation.  This    switching   action    results 

in  the  minimization  cf  power  losses  in  the  transistors,  with 
a  savings  realized  in  reduced  heat  sinking  requirements  and 
in  the  usage  of  less  expensive  power  transistors.  Since  the 
power  transistors  are  switched  on  and  off  at  a  frequency 
beyond  the  system  bandwidth,  the  motor  will  filter  the  high 
frequency  components  of  the  modulated  signal  and  respond 
only    to    the   signal's   low   frequency   components. 

Closing  a  feedback  loop  around  the  pulsewiith  modulated 
amplifier  results  in  the  amplifier  behaving  as  either  a 
current    or    a  voltage    source.  The    feedback   loop    allows    the 

motor  to  be  connected  with  the  amplifier  through  additional 
series  inductance,  resulting  in  a  smoothing  of  any  current 
ripple.  The  feedback  loap  also  allows  for  easy  current 
limiting,  merely  by  limiting  the  output  of  the  feedback 
summing  amplifier.  Finally,  a  feedback  loop  provides  output 
short  circuit  protection,  as  the  output  current  is  deter- 
mined by  the  input  voltage  without  regards  to  the  output 
impedance. 
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B.       MODULATION    TECHNIQUES 

There   are  two      basic    methods    for    obtaining      a    pulsewidth 
modulated   signal. 

1-  Ili e   Dither    Method 

The  first  technique  requires  that  the  input  signal 
(Xo)  ,  be  added  to  a  high  frequency  sawtooth  signal  (also 
known  as  a  dither  signal).  After  summing,  the  resultant 
signal  (Y1) ,  is  fed  into  a  relay  element.  The  relay  element 
converts  the  summed  signal  into  a  two  level  output  (Y2)  , 
which  then  switches  from  +V  to  -V  whenever  Y1  experiences  a 
zero    crossing,    as   shown    in    Figure   2.1. 

The  duty  cycle  (a)  of  the  output  signal  is  related 
to  the  input  signal  and  magnitude  (E)  of  the  sawtooth  signal 
by: 

a    =    (E   +    X)/2E  (eqn    2.1) 

Utilizing  this  technique  has  as  an  advantage  the  fact  that 
one  may  control  the  frequency  of  the  supplied  sawtooth 
signal  without  any  changes  to  the  motor  controller 
circuitry. 

2-  The  Limit  Cy_cle  Method 

The  second  technique  for  producing  a  pulsewidth 
modulated  signal  is  by  closing  a  feedback  loop  around  a  two 
level  switch.  The  feedback  signal  causes  the  system  to 
exhibit  high  frequency  limit  cycle  behavior.  Taking  a 
velocity  control  system  as  an  example,  if  the  velocity  error 
signal  (reference  signal  minus  the  actual  velocity)  were 
negative,  indicating  that  the  velocity  was  too  high,  the 
output  signal  from  the  switch  would  tnen  serve  to  drive  the 
velocity  lower.   As  the  velocity  of  the  system  dropped  below 
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Figure   2.1        Creation   of   a    PWM    Signal 
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the  reference  level,  the  error  signal  would  then  become 
positive,  causing  the  ouput  of  the  switch  to  change  states 
and  consequently  force  the  velocity  to  once  again  increase. 
In  this  manner  the  switch  output  would  thus  exhibit  oscilla- 
tions, describing  a  pulsetrain  whose  frequency  is  determined 
by  the  voltage  levels  at  which  tae  switch  operates  and  by 
the  dynamics  in  the  feedback  path.  The  frequency  of  the 
velocity  waveform  would  necessarily  be  the  same  as  that  of 
the  output  of  the  switch,  and  under  steady  state  conditions 
would    become   constant  and   periodic.  The    pulsewidth    of    the 

signal  at  the  output  of  the  switch  is  dependent  upon 
specific  system  dynamics.  For  tne  case  of  a  d.c.  motor 
operating  under  load,  the  pulsewidth  is  directly  related  to 
the  lead  on  the  motor.  The  specific  details  of  a  limit 
cycle   velocity    control   scheme   is    examined    in    Chapter    3. 

C.       ANALYSIS   OF    THE    PULSEWIDTH    MODULATED    SIGNAL 

A  diagram  of  an  ideal  dc  motor  is  shown  in  figure  2.2. 
The  "freewheeling"  diode  (FWD)  serves  to  bypass  the  motor 
during  the  pulse  off  period,  allowing  the  armature  current 
to   circulate.  Figure    2.3    shows      the   typical      steady    state 

current  and  voltage  relationship  in  a  pulsewidth  control 
scenario. 

Because  the  supply  voltage  is  being  switched  at  frequen- 
cies typically  on  the  order  of  magnitude  of  5  KHz,  it  is 
important  to  study  the  power  losses  within  the  motor  with 
the  power  being  pulsed  on  and  off.  A  first  approximation  to 
the  evaluation  of  power  losses  due  to  heating  within  the 
armature   resistance    is   seen    in   equation    2.2. 

?1     =   R(l2rms>  (eqn    2.2) 
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Figure   2.2        Basic    DC    Motor   Schematic 

We    may   now    define   the  current      form   factor     (k.)       as    the   ratio 
of    the   RMS    current    to  the   average   current: 


k    "    Irm£/: 


ave 


(eqn    2.3) 


Developed  motor  torgue  is  directly  proportional  to  motor 
current  (equation  2.4),  allowing  motor  losses  (under  ?KM 
conditions)    to    be   described    as   in   equation    2.5. 


Tg    -    Kt*Iave 


p.,     =    R*k2*I2 
rl         a    *      x    ave 


(eqn    2.4) 
(eqn    2.5) 


Substituting  equation  2.4  into  equation  2.5  demonstrates 
that  motor  losses  are  dependent  on  the  current  form  factor 
and    the    output    motor    torque,    as    seen    below: 


Pn    =    (R/K    2)  *k2*T    2 
It's 


(eqn    2.6) 
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Figure  2.3        PWM    Steady   State   Behavior 

The  motor  form  factor  has  a  large  influence  on  motor 
heating  and  hence,  power  losses.  Since  performance  in  speed 
control  systems  is  often  limited  by  power  dissipation 
constraints,  it  is  important  to  determine  the  form  factor 
for  a  given  PWM  scheme.  The  relationship  between  motor 
armature  losses  and  form  factor  is  shown  in  Figure  2.4 
[fief.    4], 

In  order  to  determine  the  form  factor  for  a  given 
system,  one  must  first  be  able  to  determine  a  system's 
average      and      RMS      currents.  The      differential      equations 

describing    the      motor  action   for      the    basic   dc      motor   system 
(as    seen    in    Figure    2.2)     are    as    follows: 


Pulse    on: 


Ldi/dt  =    V   -    RI    -    ku*w 


(egn    2.7) 
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Figure   2.4        Additional  Armature   Losses   vs.    Form   Factor 


Pulse   off:        Ldi/dt  =    HI    -    Kb*v 


(egn  2.8) 


where  Kt  is  the  counter  e  mf  constant  and  w  is  the  motor 
speed.  Solving  the  differential  equation  for  motor  current 
(and  referring  to  Figure  2.  3)  yields: 


Pulse  on:    1  =  1   -  (I   -  I   )exp(-t/t  )       (egn  2.9) 

n      n     oo'         n 


Pulse  off:   I  =  (lQ  ♦  lm)  exp(-t/tn)  -  I 


m 


(egn  2.  10) 
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where    tn   is   the    system  electrical    time   constant, 


I      =    (V  -   K   *w)/R  (egn    2.11) 


and 


I      =    K,  *w/R  (egn    2.  12) 

m  b 

The    average   and    RMS    currents    are   then: 


I  =    (al      -    bl    )  /E  (egn    2.  13) 

ave  v      n  m  3 


12  =    (1/p)  ((al    2  +    bl    2,-   -    t    (I      )  (I    ))  (egn    2.14) 

rms  n  m  n      od         t 


where 


a    =   t    *ln({I      -   I      )/(!      -   I    ))  (egn    2.15) 

n  n  oo  n  o 


b    =    t    *ln((I       +   I    )/(!         +   I   ))  (egn    2.16) 

n  o  m  oo  m 


Jod    =     (Jo   "   Ioo)  (e*n    2'17) 


Zt    =     (In    *   Im)  {e^n    2'18) 


and    1/P    is    the    switching   freguency. 
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As  a  means  of  demcnstra tin g  the  viability  of  ?WM  control 
of  dc  motors,  simulations  were  conducted  using  fixed  duty 
cycle  power  pulses  to  determine  the  relationship  between 
form  factor  and  motor  load  torque.  A  computer  program  was 
written  to  analyze  the  output  motor  current  waveshape  for 
the  average  and  rms  currents  utilizing  relationships 
detailed  in  this  section.  The  program  used  to  compute  these 
currents  may  be  found  in  Appendix  C.  Figures  2.5  and  2.6 
show  that  the  form  factor  rapidly  approaches  unity  as  the 
load  on  the  motor  is  increased,  indicating  that  the  motor  is 
experiencing  only  slight  additional  losses  (in  terms  of 
percentages)  due  to  the  power  pulsing  effect  as  compared 
with    a      constant   voltage     supply   arrangement.  It   is      also 

pointed  out  that  the  motor  runs  more  efficiently  (lower  form 
factor)  at  higher  frequencies  for  a  given  load  torque.  This 
is  due  to  the  fact  that  as  the  motor  is  pulsed  more 
frequently,  the  motor  speed  will  not  drop  off  as  fast  and 
hence,  the  energy  required  to  move  the  mass  of  the  rotor 
back  to  its  steady  state  speed  will  not  be  as  great. 
However,  switching  losses  in  the  transistors  will  usually 
limit    the    switching    frequency    to    less    than    ten   KHz. 
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FORM  FACTOR  VS.  LOAD  TORQUE 
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Figure  2.5        Form    Factor    vs.    Load    Torque    (1    KHz) 
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FORM  FACTOR  VS.  LOAD  TORQUE 
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III.     PULSEWIDTH    MODULATED    SPEED    CONTROL 

The  ability  of  a  dc  motor  to  maintain  a  given  speed  when 
a  load  torque  is  applied  is  generally  referred  to  as  speed 
regulation.  Although  a  dc  motor  by  itself  is  an  open  loop 
system,  the  prescence  of  the  back  electromotive  force  (temf) 
signal  serves  to  close  a  natural,  "built-in",  feedback  loop, 
as  shown  in  Figure  3.1  However,  because  the  dc  motor  is 
intrinsically  an  open  loop  system  with  relatively  constant 
power  input,  as  the  load  torgue  increases,  the  speed  will 
decrease,  and  hence,  no  speed  regulation  may  be  achieved. 
In  order  to  maintain  a  constant  speed,  the  input  power  must 
increase   with   the  applied   load. 
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Figure   3.1         Block   Diagram   of   a   Basic    DC    Motor 

Figure      3.2      illustrates      the        basic      principle      behind 
pulsewidth    modulated    amplifiers.         flhat    is    important    to   note 
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is  that  as  the  load  en  the  system  increases,  the  duty  cycle 
or  pulsewidth  of  the  input  signal  also  increases.  One 
thought  then,  is  to  attempt  motor  speed  control  by  deter- 
mining wiiich  motor  parameters  are  changing  relative  to 
varying  loads  and  to  then  make  pulsewidth  a  function  of  one 
or  more  of  those  parameters.  Because  motor  current  varies 
linearly  with  load  torgua,  some  form  of  current  feedback 
appears  to  be  the  logical  selection  for  implementation  of  a 
speed  control  scheme.  A  specific  current  feedback  control 
technique  was  investigated,  as  well  as  a  limit  cycle  control 
method,  the  details  of  which  will  now  be  presented. 
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Figure  3.2    Pulse  Width  as  a  Function  of  Load  Torgue 
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A-       COEHZNT    FEEDBACK 

Prior  to  developing  any  specific  control  schemes, 
studies  were  conducted  of  motor  behavior  relative  to  a  fixed 
pulsewidth  power  pulsed  input,  at  a  frequency  of  5  KHz. 
Motor  speed  was  found,  as  in  the  case  of  constant  supply 
voltage,  to  vary  linearly  with  load  torque  ,  as  shown  in 
Figure  3.3.  As  an  aside,  all  studies  were  conducted  with 
the  load  torques  ranging  from  zero  to  eighty  ounce-inches, 
as  this  range  represented  the  linear  range  of  operation  for 
the  motor  modelled  in  the  Thomas  study  [Ref.  2].  The  plot 
of  average  current  vs.  load  torque  for  the  fixed  pulsewidth 
simulations   are    shown   in   Figure   3.4. 

Studying  the  curves  found  in  Figures  3.3  and  3.4  led  to 
the  conclusion  that  a  scheme  for  speed  control  of  the  motor 
could  be  developed  with  the  pulsewidth  being  made  directly 
proportional  to  the  average  motor  current.  The  basic  form 
of   the   motor  pulsewidth   was    decided   to    be   as    follows: 

PW    =      DCF   +    Tave    *    K  (egn    3.1) 

The    motor   parameters   found    in   equation   3.  1    may   be    defined   as 
follows : 

PW  =   the  input    pulsewidth    {duty    cycle) 

Iave  =   average   motor  current 

DCF  =   a    dc   term    which   establishes    no    load   speed 

K  =   the  current    feedback   constant 

It  was  decided  that  the  motor  would  run  at  a  minimum  of 
50%  duty  cycle  pulses  to  minimize  power  losses  which  would 
occur  at  smaller  duty  cycles  in  light  load  conditions.  With 
a  frequency  set  at  5  KHz,  this  necessarily  fixed  no  load 
speed  at  approximately  1375  rpm.  Since  large  current  tran- 
sients  could  be    expected   when    the      motor   was    started    or   when 
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AVERAGE  SPEED  VS.  LOAD  TORQUE 


o 

© 
o 

C3 


o 

© 

CM 


□ 

50% 

LEGEND 
PULSE  DUTY 

CYCLE 

0 

75% 

PULSE  DUTY  CYCLE 

A 

NO  PWM 

0.0        10.0        20.0        30.0        40.0        50.0        60.0        70.0        80. 


LOAD  TORQUE  (OZ-IN) 


Figure   3.3        Motor    Speed   vs.    Load   Torque    (fixed    pulse    width) 
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MOTOR  CURRENT  VS.  LOAD  TORQUE 
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Figure   3.U        Current   vs.     Load   Torque    (fixed   pulse    width) 
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the  load  condition  was  changed,  a  limiter  was  written  into 
the  simulation  program  to  set  the  maximum  pulsewidth  at  just 
less  than  the  100%  duty  cycle  point.  It  was  not  set  at  100% 
due  to  difficulties  encountered  with  the  simulation 
language.  Other  additions  to  the  basic  computer  model 
included  a  limiter  to  prevent  pulsewidths  with  less  tnan  505 
duty  cycles  as  well  as  a  current  limiter  to  prevent  negative 
currents.  The  latter  was  added  to  simulate  the  effect  of 
the  addition  of  the  "freewheeling"  diode  described  in  the 
preceding  chapter. 

To  set  the  no  load  speed  at  1375,  it  was  necessary  to 
obtain  parameters  that  would  establish  PW  equal  to  .5  at 
zero  load  torque.  The  value  of  K  was  determined  from 
Figures  3.3  and  3.4  by  noting  that  the  75%  duty  cycle  condi- 
tion, at  approximately  1350  rpm,  occurred  at  a  load  torque 
of  40.0  oz-ins,  which  also  corresponded  to  an  average  motor 
current  of  2.50  amperes.  Setting  PR  to  .75  and  DCF  to  .48 
in  egn.  3.1  led  to  a  K  value  of  .107.  A  value  of  DCF  of  .48 
results  in  an  approximately  50%  pulsewidth  modulated  signal, 
as  no  load  motor  current  is  approximately  .002  amps.  It  was 
felt  that  higher  speeds  could  then  be  achieved  by  increasing 
DCF,  as  the  basic  relationship  between  pulsewidth  and  speed 
appeared   linear  for  any  given  load  torque. 

Extensive  simulations  were  conducted  utilizing  the  feed- 
back control  relationship  as  shown  in  equation  3.2.  Graphic 
results  for  these  simulation  trials  are  shown  in  Figure  3.5. 

Pfl  =  .48  +  (I    *.  107)  (egn  3.2) 

v  ave  \    i  * 

It  is  clearly  evident  from  Figure  3.5  that  the  control 
scheme  utilized  performed  u nsatisfactoriliy  for  its  task  of 
maintaining  constant  speed  throughout  the  given  range  of 
load  torques.  While  certain  variations  might  have  been 
expected     in  the      output      speed,         the    results      demonstrated 
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nonlinearities  in  motor  performance  wnich  were  clearly 
unsuitable  for  its  given  application.  One  of  the  major 
problems  encountered  with  this  control  scheme  stemmed  from 
the  fact  that  the  average  current  was  used  as  a  feedback 
parameter,  rather  than  the  actual  motor  current (the  ripple 
present  in  the  motor  current  was  deemed  to  be  too  high  to  be 
used  in  a  velocity  control  scheme  based  on  current  feed- 
back) .  The  average  motor  current  proved  to  be  unsatisfac- 
tory for  the  given  task  for  two  reasons:  1)  average  motor 
current  requires  time  to  approach  the  system's  actual 
average  current  value  due  to  the  changes  to  the  current 
incurred  by  transients  such  as  are  caused  by  changes  to  the 
system's  steady  state  behavior  and  thus  adds  significantly 
to  the  motor's  settling  time  to  variances  in  load  or 
commanded  speed,  and  2)  apparent  nonlinearities  which  would 
appear  if  pulse  duty  cycle  were  plotted  as  a  function  of 
load  torque,  which  runs  contrary  to  the  initial  assumptions 
upon    which    the    control   scheme    was    devised. 

B.       VELOCITY  LIHII    TECHHIQOE 

The  velocity  limit  control  scheme  has  been  developed 
under  the  assumption  that  a  pure  velocity  command  has  been 
issued  by  the  motor  control  logic.  It  is  recognized  that 
other  systems  might  issue  torgue  commands  to  the  motor  in 
response  to  a  generated  missile  fin  position  error  signal 
and    the      current   motor  speed.  Simulation   of      the    complete 

electromechanical  actuator  is  not  the  intent  of  this  thesis, 
and  hence,  will  not  be  attempted  here.  The  schematic 
diagram  for  the  network  used  to  implement  the  limit  control 
scheme  is  shown  in  Figure  3.6.  This  network  will  control 
the  motor  speed  in  such  a  way  that  if  speed  is  below  its 
commanded  level  (plus  a  pre-defined  tolerance),  the  power  to 
the    motor   is  turned    on.      If    the    motor    speed   rises    above    this 
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Figure   3.5        Motor    Speed   vs.    Load   Torque    (current    feedback) 
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set    point,    the    power    is   turned   off.       The    controlled    velocity 
waveshape    is  diagrammed    in    Figure    3.7. 
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Figure   3.6         Schematic   Diagram   Of    Pulsewidth    Modulator 

Operation  of  the  system  in  the  limit  mode  is  relatively 
straightforward.  Referring  to  Figure  3.7,  if  the  motor 
speed  is  below  (VCOM  +  VTOL) ,  transistor  QM  is  switched  on, 
allowing  motor  current  to  flow.  When  the  speed  reaches 
(VCOM  +  VTOL)  ,  QM  is  switched  off,  which  then  induces  a 
large  voltage  across  the  inductor  terminals,  owing  to  a 
rapid  rate  of  change  in  the  inductor  current.  This  induced 
voltage  turns  on  diode  DM,  which  provides  a  path  for  the 
decaying  motor  current.  When  motor  speed  decays  past  (VCOM 
-  VTC1),  QM  is  switched  back  on  again.  Transistor  QB  and 
diode  DB  are  utilized  when  the  motor  is  operated  in  the 
regenerative  braking  mode  of  operation.  In  this  mode,  the 
dc    motor   is    used   as    a   generator      as    the    inertia    of    the    rotor 


31 


MOTOR 
SPEED 


VCOM+  VTOL 


VCOM, 


VCOM-VTOL 


time 


Figure   3-7        Motor    Velocity   Waveshape 

is  converted  into  electrical  energy  to  charge  the  system's 
dc  voltage  supply  while  controlling  the  braking  of  the 
rotating   masses.  Again,       how    this      mode    is    utilized      is   a 

function  of  the  design  of  the  motor's  electromechanical 
actuator,      and       will    not      be   further      commented    upon.  The 

control  of  transistors  QB  and  QH  can  be  accomplished 
utilizing  voltage  comparators,  where  one  input  is  a  voltage 
proportional  to  the  commanded  speed  plus  the  allowed  speed 
tolerance,  while  the  other  input  is  a  voltage  proportional 
to      the    actual      motor     speed.  What    becomes      important      to 

realize  now  is  that  the  pulsewidth  and  the  pulse  frequency 
are  both  variable,  and  will  be  dependent  on  certain  dynamics 
of    the   system. 
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A  trial  simulation  was  conducted  using  the  limit  control 
scheme.  The  speed  tolerance  was  set  at  five  rpm,  and  the 
commanded  speed  was  set  to  1400  rpm.  Table  I  contains  the 
results  of  the  simulation.  This  data  is  also  represented 
graphically  in  Figure  3.8.  Loading  of  the  motor  was  accom- 
plished using  a  terminated  ramp  signal;  the  terminal  value 
of  the  ramp  is  the  desired  motor  loading.  The  speed  accu- 
racy of  the  motor  is  defined  as  the  difference  between  the 
maximum  and  minimum  motor  speed  divided  by  the  commanded 
speed. 

TABLE   I 
Hotor  Characteristics 


Load      Torgue  Speed  (rpm)  %   accuracy 

(oz-in)  maximum  minimum 

0.0 

16.0 

32.0 

48.0 

64.0 

80.0 


The  results  of  the  initial  simulations  indicates  that 
positive  control  of  the  commanded  motor  speed  may  be 
accomplished  utilizing  the  limit  cycle  method.  All  further 
studies  are  therefore  based  on  a  computer  model  whose  speed 
is   controlled   in    this   manner. 

Studies  of  system  characteristics  (pulsewidth,  pulse 
frequency,  speed  regulation)  are  presented  in  detail  in  the 
next  chapter  in  order  to  better  define  the  operational 
enevelope  of  the  modelled  dc  motor  using  the  velocity  limit 
control  scheme.  Addittionally,  the  effects  of  the  addition 
of  series  inductance  are  also  investigated  to  determine  its 
affects   in    regards    to   ripple   reduction. 
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MOTOR  SPEED  ACCURACY 
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17.     PERIOMANCE    OPTIMIZATION 

The  simulations  described  at  the  end  of  the  previous 
chapter  were  not  necessarily  indicative  of  the  optimum 
performance  potential  of  the  modelled  dc  motor.  It  appears 
evident  that  there  must  exist  means  by  which  the  current 
ripple  and  the  speed  regulation  may  be  respectively  reduced 
and  improved.  One  way  to  achieve  improved  speed  accuracy 
(which  implies  that  the  ripple  present  in  the  motor  velocity 
is  reduced)  is  to  decrease  the  speed  tolerances  established 
for    motor   operation.  The    effects   of    the      variation    of   the 

speed  tclearnce  settings  will  be  examined  shortly.  However, 
to  reduce  the  current  ripple,  which  in  turn  translates  to 
reduced  power  losses  in  the  motor,  it  has  been  suggested 
that  one  must  add  inductance  in  series  with  the  motor  arma- 
ture [fief.  5].  The  effect  of  the  addition  of  series  induc- 
tance on  motor  performance  is  studied  in  the  following 
section. 

A.       ADDITION   OF    SERIES   INDUCTANCE 

A  dimensionless  current  ripple  may  be  defined  as  the 
curent  ripple  multiplied  by  the  motor  armature  resistance 
and      divided  by      the    supply      voltage.  The   current      ripple 

itself  is  defined  as  the  difference  between  the  current  at 
the  time  the  motor  is  pulsed  on  and  the  time  when  the  motor 
is  pulsed  off,  or  simply  as  the  difference  between  the 
minimum   and   maximum    currents,    as    follows: 

I    =    i(aT)    -    i(0)  (eqn    4.1) 
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Foe  th€  case  where  only  one  power  supply  is  used  (unidirec- 
tional drive),  dimensionless  current  ripple  (Id)  may  be 
defined    as    follows; 

Id    =    (1  +  exp(-tn)-exp(-a*tn)-exp((1-ajtn))/g  (eqn    4.2) 

where: 


tn   =   L/E  (eqn    4.3) 


and 


g    =    1    -   exp(-tn)  (eqn   4.4) 

Dimensionless  current  ripple  may  be  plotted  versus  the 
puisewidth,  or  duty  cycle,  with  the  ratio  of  the  period  of 
the  PWM  signal  to  the  electrical  time  constant,  tau,  to  form 
a  family   of   curves,    as   in  Figure    4.1. 

As  can  be  seen  from  Figure  4.1,  the  magnitude  of  the 
current  ripple  depends  to  a  great  extent  upon  the  ratio  of 
the  pulse  period  to  the  motor  electrical  time  constant.  We 
should  therefore  expect  that  that  reduction  of  the  motor 
current  ripple  can  be  accomplished  through  the  addition  of 
series   inductance,    which   reduces    the   magintude   of   tau. 

The  effect  of  adding  series  inductance  to  the  motor  may 
be  seen  by  writing  the  differential  equations  for  the 
simplified  circuit  diagram  of  a  dc  motor  as  shown  in  Figure 
4.2.  Applying  Kirchoff's  Law  and  summing  voltages  around 
the  lcop  results  in  a  cucr ent-voltage  relationship  as  shown 
in    equation   4.5: 

V    =     (1   +   L1)di/dt    +   Ri   +    Kb-w  (eqn    4.5) 

Taking  the  Laplace  transform  of  equation  4.5  yields: 
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Figure  4-1    Current  Ripple  vs.  Pulse  Duty  Cycle 
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V  (s)     =    s  (L    +  L1)I(s)     +    RI  (s)     +    Kh-W(s) 


(eqn    4.6) 


Solving    equation    4.6    for   motor   current    yields: 


I(s)    =    (l/s(L    +    LI)    +    R)).(V(s)    -    KbW(s)) 


(eqn    4.7) 


From    equation   4.7   the  electrical    time    constant   is: 


tn    =     (L    +    L1)/R 


(eqn    4.8) 


Figure  4.2    Basic  DC  Motor  Circuit  Diagram 

It  is  tc  be  noted  that  the  previous  derivation  ignores 
the  affect  of  the  resistance  which  will  accompany  the 
additional  series  inductance.  However,  it  will  be  assumed 
that  a  value  for  the  electrical  time  constant  may  be  set  as 
in  equation  4.8  through  careful  selection  of  the  type  and 
quantity  of  series  inductance  added  to  the  motor. 
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To  test  the  effects  of  the  addition  of  series  induc- 
tance, the  computer  model  was  modified  to  reflect  the  addi- 
tion of  inductance  egual  in  magnitude  to  the  inductance 
present  in  the  motor's  windings  (ignoring  the  change  in  the 
total  resistance  for  the  reasons  mentioned  above) . 
Simulations  were  then  conducted  throughout  the  load  torgue 
range  of  from  zerc  to  eighty  ounce-inches  with  motor 
velocity  set  at  at  1400  rpm  (approximately  half  of  the 
motor's  no  load  speed)  and  speed  tolerance  set  to  ±1  rpm. 
Table  IT  summarizes  the  results  of  these  simulation  trials. 
The  data  for  simulations  made  without  the  additional  induc- 
tance is  shown  in  Table  III  and  is  included  for  sake  of 
comparison    with    the    data   in    Table   II. 

TABLE   II 
Inductance  Effects   on   Motor  Operation 


Load    Tor 

gue 

lave 

Irms 

Form    Factor 

Current 

(oz-ins) 

(amps) 

(amps) 

Ripple    {%) 

05.0 

0.2  57 

0.454 

1.  76  5 

.1100 

16.0 

0.9  95 

1.  131 

1.  138 

.0874 

32.0 

1.967 

2.0  32 

1.033 

.0948 

48.0 

3.0  23 

3.064 

1.012 

.0562 

64.0 

4.009 

4.028 

1.005 

.0787 

80.0 

5.011 

5.031 

1.  004 

.  1295 

TABLE    III 
Motor   Performance     (no   series   inductance) 


Load    Torgue    xave 
(oz-ins)  "  (amps) 


05.0 

16.0 
32.0 
48.0 
64.  0 
80.0 


0.  179 
1.019 

1.  983 
3.  013 
4.028 
5.017 


Irms 

(amps) 

0. 

.347 

1. 

.10  1 

2. 

.026 

3. 

.370 

4. 

.040 

5. 

.026 

Form  Factor 


1.939 
1.081 
1.022 

1.008 
1.003 
1.002 


Current 
Ripple  (%) 

.326 
.259 
.279 
.261 
.226 
.285 
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Figure   4. 3        Inductance   Effects    on   Current    Ripple 
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Examination  of  the  data  in  Tables  II  and  III  demon- 
strates clearly  the  ripple  reduction  obtained  through  the 
addition  of  the  series  inductance.  Motor  form  factor  was 
calculated  to  show  that  pwer  losses  (at  low  load)  in  the 
transistors  are  also  reduced,  and  tails  represents  an  addi- 
tional benefit  gained  from  the  added  inductance.  However, 
while  series  inductance  will  reduce  the  current  ripple 
without  affecting  steady  state  behavior,  it  will  have  an 
affect  on  the  transient  behavior  of  the  system,  as  will  now 
be  shown. 

Assuming  now  that  the  system  is  operating  in  steady 
state,  a  step  input  command  (such  as  a  change  in  the  motor's 
commanded  velocity)  will  force  the  pulsewidth  modulated 
signal  to  the  full  on  condition.  The  response  of  the  system 
will  now  be  limited  by  the  motor's  electrical  time  constant. 
In  the  case  where  series  inductance  has  been  added  to  the 
system,  the  response  of  the  motor  to  a  step  input  will  be 
slower  than  if  the  additional  series  inductance  were  not 
present. 

To  test  this  effect,  simulation  trials  were  conducted 
where  the  motor  was  allowed  to  achieve  a  steady  state  speed 
of  1000  rpm  and  then  was  subjected  a  step  input  command  to 
increase  motor  speed  to  140  0  rpm.  Toruue  load  was  32  oz-ins 
and  the  speed  tolerance  setting  was  ±  1.  The  response  time, 
or  the  time  required  for  the  motor  to  settle  at  the  new 
commanded  speed,  was  measured  for  trials  in  the  motor's 
standard  conf iguraticn  and  for  the  case  where  series  induc- 
tance was  added.  For  the  first  case,  where  there  was  no 
additional  inductance,  the  response  time  was  measured  at 
approximately  2.5  milliseconds.  When  series  inductance  was 
added  to  the  system,  the  response  time  slowed  to  approxi- 
mately five  milliseconds.  Thus  if  one  decides  to  add  series 
inductance  to  the  motor  to  reduce  tne  current  ripple 
effects,  that  decision  must  be  tempered  by  the  fact  that  the 
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transient  response  cf  the  system  will  change  due  to  the 
change   in   the   motor's  electrical   time   constant. 

B.       REDUCTION    OF    VELOCITY    HIPPIE 

Optimization  of  motor  performance  will  reguire  that  the 
ripple  content  of  the  motor  velocity  at  steady  state  tie 
minimized.  To  reduce  the  ripple,  it  will  then  be  necessary 
to  reduce  the  speed  tolerance  settings  which  will  establish 
the    pulsewidth       modulated   input      signal    to      the    motor.  To 

determine  the  effects,  if  any,  on  motor  performance  (other 
than  the  reduction  of  velocity  ripple),  it  was  necessary  to 
perform  a  number  of  simulations  with  varied  speed  tolerance 
settings. 

Three  simulation  trials  were  conducted,  with  the  speed 
tolerance   set   at      ±    5  rpm,      ±1      rpm   and    ±. 1    rpm.  Table   IV 

summarizes  the  results  of  the  simulations.  The  data  for 
Table    IV   is   plotted    in   Figures   4.4   and   4.5. 

TABLE   I? 
Performance  Trials    for  Various   Speed   Tolerance    Settings 


Speed  Tolerance    (rpm) 
±5  ±1  ±.1 

Load    Torque     ripple      form  ripple      form  ripple 

(oz-in) "  (%)  factor         (>)  factor  (%] 

05.0 
16.0 
32.0 
48.0 
64.0 
80.0 


It  is  apparent  that  reducing  the  speed  tolerance  does 
indeed  reduce  the  ripple  content  of  the  motor's  speed. 
Additionally,  as  the  tolerance  is  reduced,  so  does  the  fore. 
factor  of  the  motor(for  specific  loading),  indicating  higher 
motor    performance      efficiency.         Thus   as      speculated    before, 
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2.  19 

2.860 

0.89 

1.765 

.214 

3.00 

1.819 

1.16 

1.  138 

.242 

3.44 

1.372 

1.22 

1.  033 

.235 

3.57 

1.283 

1.21 

1.012 

.228 

3.71 

1.095 

1.26 

1.005 

.228 

3.79 

1.054 

1.29 

1.003 

.250 

f 

orm 

factor 

1 

.7  09 

1 

.121 

1 

.026 

1 

.011 

1 

.0  05 

1 

.002 

with  lower  speed  tolerances,  tnere  exists  improved  motor 
performance.  Of  course  this  should  have  been  intuitively 
obvious  even  prior  tc  conducting  the  simulation  trials  due 
to  the  fact  that  the  greater  the  speed  tolerances,  the 
greater  the  inertia  obtained  by  the  rotating  mass  of  the 
rotor  before  a  speed  limit  is  reached  and  thus  greater 
amounts  of  energy  would  have  to  be  expended  in  order  to 
increase  the  motor's  speed  back  to  the  current  commanded 
speed. 

C.       CCMHENTS 

We  have  seen  where  the  simulated  performance  of  a 
modelled  brushless  dc  motor  has  been  improved  through  the 
addition  of  series  inductance  and  the  optimizing  of  the 
speed  tolerance  settings  of  the  chopper  control  logic. 
There  are  reasons  to  assume  that  the  performance  enhance- 
ments noted  in  this  chapter  may  not  necessarily  be  realized 
in      an      operational      cruise    missile      scenario.  The      final 

chapter  of  this  report  discusses  where  caution  need  be  taken 
when  reviewing  this  work  and  prior  to  applying  these  results 
in    the   use    of   brushless   dc    motors   in   aerospace   applications. 
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VELOCITY  RIPPLE  STUDIES 
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Figure   4. 4        Motor   Velocity    Ripple    vs.    Load    Torque 
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V-  RECOMMENDATIONS  FOR  FURTHER  STUDIES 

As  stated  in  Chapter  3,  the  assumption  was  made  that  the 
missile's  control  logic  would  generate  velocity  commands  to 
the  motor  controller  unit.  This  may  not  necessarily  be  true 
for  all  missile  applications.  Certain  controllers  may 
generate  torque  commands  (or  equivalently,  current  commands, 
as  current  is  directly  proportional  to  motor  torque)  to  the 
motor  to  provide  a  constant  torque  to  affect  the  required 
missile  maneuvers.  The  effect  that  this  would  have  on  the 
simple  control  scheme  studied  in  this  report  could  prove  to 
be  fairly  significant,  and  will  now  be  looked  at  in  closer 
detail. 

A.   BEYOND  SPEED  CONTROL 

One  major  change  that  would  be  required  of  the  motor 
controller  if  torque  commands  are  to  be  issued  from  the 
missile  control  logic  is  that  both  position  and  current 
loops  would  necessarily  have  to  be  closed  around  the  PWM 
amplifier.  Of  course,  the  implication  of  generating  torque 
commands  is  that  the  operating  envelope  of  the  motor  by 
necessity  would  have  to  be  defined  to  include  the  plugging, 
braking  and  the  regenerative  braking  modes  of  operation. 
The  importance  inherent  within  the  inclusion  and  modelling 
of  these  modes  is  best  justified  by  noting  that  the  regener- 
ative braking  mode  serves  not  just  to  control  the  braking 
torque  applied  to  the  rotor  but  also  to  recharge  the 
missile's  supply  batteries.  The  model  would  then  need  to 
have  logic  blocks  which  could  recognize  when  each  operating 
mode  would  be  appropriate  and  then  generate  the  requisite 
commands  to   the  motor  control   logic.    The   idealized  step 
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response  of  the  fin  actuator  system  to  command  chance  in  fin 
position  is  shown  in  Figure  5.1  with  the  required  motor 
operating    modes    shown  as   a    function    of    flap   position. 

Since  1)  a  position  loop  has  yet  to  be  closed  and  2)  no 
accounting  has  been  made  for  the  modes  of  operation  ether 
than  motoring  in  the  forward  direction/  it  is  most  highly 
recommended  that  these  areas  be  investigated  to  more  accu- 
rately model  the  motor  as  an  integral  part  of  an  overall 
missile  fin  actuation  system.  It  is  emphasized  here  that 
modelling  the  position  control  of  the  fin-motor  system  would 
represent  the  next  logical  step  in  accurately  simulating  the 
dynamics   of   the    electromechanical   actuator    sytem. 

B.       SIMULATION    DEFICIENCIES 

Studies  of  the  modelled  motor  have  been  concentrated  on 
the  response  of  the  system  to  step  input  commands  under 
constant  load  conditions.  Of  course,  under  normal  operating 
conditions,  the  system  most  undoubtably  would  be  subject  to 
a  series  of  varied  load  conditions  as  the  missile  steers  its 
course  towards  the  target.  The  load  on  the  motor  would  then 
be  a  function  of  the  aerodynamics  to  which  the  missile  is 
subjected  during  flight,  i.e.,  missile  speed,  attitude, 
acceleration,  etc.  In  order  then  to  model  the  system  better 
under  the  dynamics  of  flight,  the  effort  to  close  the  posi- 
tion loop  should  be  followed  by  more  precise  load  studies, 
so  that  the  motor's  behavior  may  be  studied  within  a  context 
more  closely  related  to  its  predicted  operating  environment. 
Other  studies  that  would  prove  worthwhile  include  (but 
certainly  are      not    limited    to)  :  the   study   of      the   voltage 

switching  affects  on  the  controller's  power  transistors, 
designing  the  hardware  required  to  implement  the  notor 
controller  logic,  as  well  as  the  associated  software. 
Studying   of    the    effects   of    closing    a    phase-locked    servo   loop 


47 


Ffn  Angle 
Regenerative  Braking  Mode 


Motoring  (-)  Mode 


Motoring  (*)    Mode 


Comanded  Flap 
Position 


Time 


Figure   5. 1        Theoretical   Fin   Actuator   Response 
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for    velocity     control   might      also   be      considered    due      to    the 
system's    precise  speed   regulation   capaoility. 

One  caution  worth  noting  is  that  the  model  used  as  a 
basis  for  all  simulation  studies  assumed  a  linear,  average 
flux  back  emf  waveform.  Ihe  back  emf  signal  is  in  actuality 
siuusoidal  in  nature,  and  would  therefore  result  in  the 
addition  of  fundamental  and  harmonic  frequency  components  to 
the  motor  parameters  which  were  studied,  sucn  as  current  and 
velocity.  It  is  recommended  that  for  further  studies  a  more 
advanced  motor  model  which  simulates  sinusoidal  back  emf  be 
utilized. 

C.       SUMMARY   OF    RESULTS 

Pulsewidth  modulation  has  been  shown  to  be  a  viable 
method  of  accurately  and  reliably  controlling  the  velocity 
of  a  brushless  dc  motor.  Form  factor  studies  indicated  that 
power  losses  in  the  switching  transistors  may  be  minimized, 
thus  allowing  for  smaller  power  transistors  and  reduced  heat 
sinking,  which  translates  into  reduced  costs  in  controller 
design  and  implementation.  The  addition  of  series  induc- 
tance was  shown  to  a  have  a  very  definite  impact  on  motor 
current  ripple,  reducing  it  significantly  in  comparison  to 
simulations  conducted  without  the  series  inductance.  Use  of 
the  limit  cycle  pulsewidth  modulation  scheme  was  shown  to  be 
a  superior  method  for  implementing  pulsewidth  modulation. 
Areas  where  further  research  efforts  may  continue  have  been 
presented,  including  specific  recommendations  for  follow-on 
studies    to   this    thesis. 
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APPENDIX    A 
THE    MOTOR    MODEL 

The  motor  that  was  modelled  was  a  commercially  available 
brushless   dc  motor.  The    current    and   speed      carves    for    the 

motor  have  been  shown  in  Figures  3.3  and  3.4.  The  motor  is 
a  three  phase,  four  pole  machine  with  the  commutation  being 
accomplished  electronically  utilizing  a  set  of  three  Hall 
effect  position  sensor  devices  and  a  set  of  six  switching 
transistors.  The  switching  takes  place  every  30  degrees  of 
mechanical  rotation. 

The  back  EMF  signal  was  assumed  to  be  directly  propor- 
tional to  the  motor  speed.  The  actual,  sinusoidal  nature  of 
the  waveshape  was  not  taken  into  account  in  the  model  used. 
What  follows  is  a  brief  description  of  certain  procedures 
that  were  added  to  the  basic  motor  simulation  program,  a 
description  of  program  variables  added  to  the  initial 
program  as  well  as  observations  made  concerning  the  execu- 
tion   of   the   simulation   program. 

A.       PBOCEDOEES    ADDED    TO    THE    BASIC    PEOGRAM 

1 .      Procedure  I  CLIP 

Procedure  I CLIP  was  added  to  the  model  to  account 
for  the  addition  of  the  freewheeling  diode  FWD.  Currents 
were  then  necessarily  clipped  (negatively  going)  at  zero 
amperes,  thus  preventing  the  circulation  of  negative 
currents. 

2  .      Procedure   VCIP 

This  procedure  takes  as  an  input  the  velocity  error 
signal      and  yields      the      motor's      input    voltage      based      upon 
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commanded  speed  and  the  established  speed  tolerance.  It  is 
within  this  procedure  that  the  limit  cycle  behavior  of  the 
system  is  established. 

3.   Procedure  RESET 

As  the  output  of  the  integrator  block  which  yields 
rotor  position  counts  upward  continuously  from  zero,  a 
procedure  was  required  to  reset  the  rotor  position  to  zero 
degrees  when  360  degrees  of  rotation  of  tne  rotor  was 
achieved.  F.ESET  keeps  track  of  the  number  of  times  that  the 
rotor  turns  past  the  360  degree  point  and  uses  this  informa- 
tion to  update  the  variable  THEST,  which  cycles  ranges  from 
zero  to  360  degrees.  The  importance  of  this  procedure  will 
ultimately  be  realized  when  a  position  loop  is  closed  within 
the  system. 

B.   MOTOR  PARAMETERS 

The  parameters  which  follow  are  those  added  to  the  basic 
program  to  achieve  the  required  speed  control  effects. 

1.  VCMD 

VCMD  is    the    commanded   motor   velocity. 

2.  VTOL 

VTOl  is  the  velocity  tolerance  limit  in  rpm. 

3.  ICLIP 

ICLIP  is  the  motor  current  which  has  been  adjusted 
to  prevent  negative  current  flow. 

4.  THE  ST 

THRST  is  the  rotor  position  in  degrees  which  ranges 
from  zero  to  360  degrees. 
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C.        NOTES    ON   PROGRAM    EXECUTION 

A  fixed  interval  integration  technique  was  used  in  lieu 
of  the  variable  step  Runge-Kutta  method  supplied  as  a 
default  integration  technique  in  CSMP-III.  The  trapezoidal 
technique  was  used  as  it  demonstrated  itself  to  perform  as 
accurately  as  the  varuable  step  methods  but  used  better  than 
50%  less  computer  time.  The  integration  interval  was  chosen 
as    .000001    seconds. 

Since  the  switching  frequency  of  the  motor  was  on  the 
order  of  2500  KHZ,  to  accurately  observe  motor  behavior 
during  the  pulse  on  and  off  periods  a  print  interval  of  less 
than  50  microseconds  (typically  20  microseconds  was  chosen) 
was  required.  Because  of  this,  it  was  difficult  to  observe 
the  microscopic  detail  of  motor  operation  in  terms  of  the 
variances  present  in  the  current  and  speed  for  periods  of 
greater  than  one  second,  as  the  CSMP  program  is  limited  to 
approximately  5500  lines  of  output.  For  the  studies  made 
for  this  report  the  limitation  encountered  did  not  pose  a 
major  problem,  but  could  prevent  an  obstacle  to  further 
studies.  Of  course,  as  the  studies  of  this  system  advances, 
the  requirements  for  such  detailed  assessment  of  motor  oper- 
ation may  not  be  present,  and  larger  print  intervals  may  be 
used,    thus    allowing    studies    to   be    of    greater   duration. 
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APPENDIX    B 
CSMP    SIMULATION    PROGRAM 


//ASKCSMP2     JOB     (1102,0116)  , ' ASKINASAAA A' ,CLASS  =  C 

//♦MAIN    OR  G=N  P  GViMl.  110  2P 

//    EXEC    CSCFXV 

//X.COKPRINT    DC    GbtAM 

//X.SYSPRINT    DC    CUMA-Y 

//X.SYS1N    DO    * 

INITIAL 

CONSTANT    KT   =     15.9.    BM    =    0.00015,     BL    =    0.0,    JL    =   O.C, 

N    =     1.0,    JM   =     C.001,     KB   =    0.112,    TLP    =    104. 
PARAMETER    LA    =    .3016,     RA   =    2.740.     ... 

TLF    =    05.,    VTCL    =    1.  i     VCMC   =    1400. 

*  Al    =    LA/PA      --       ThE    INVERSE    ELECTRICAL    TIKE    CCNlST/NT 

*  A2    =    J/B       —       THE    INVERSE    MECHANICAL    TIME    CONSTANT 

*  VTOL    —    ThE     VELOCITY     LIMIT    SETTING 

*  VCMD    —    THE     COMMANDED    MOTCR    VELOCITY 

*  TLF    —   LOAD    TCRCUE    IN    CZ-INS 

*  KB    —    EACK    EM    CCNSTANT 

*  KT    —    TuRCUE    CCNSTANT 

*  N    —    GEAR    RATIO 

*  JM    —    ROTCR    *CJ«EN7   OF    INERTIA 

*  TLP    —   LCAD    TCRCUE    WHERE    PEAK    POWER    COT    CCCJPS 

NO  SORT 

BLP    =    BL/(N**iJ 

JLP    =    JL/<N**2) 

J    =    JM   +     JLP 

B    =    £M    +    ELP 

Al    =    2.0    *    LA    /    RA 

A2    =    J    /    E 

THRST    =    0.0 

JFAC    =0.0 

DYNAMIC 

VERR    =   VC*D    -    kMRFM 

VIF=    VCLP   *    S1EPIC.0J 

TL1    =    RAMP(O.C) 

TL2    =    RA^P(.Cl) 

TL    =    TLF    *    ( TL1    -   TL2)/. CI 

VIB    =    0.0 

VI N    =    VI F    -*■     VIE 

VIN1    =    V  IN    -    VEMF 

VIN2    =    VIN1    *     (1.0/RA) 

IM    =    REALPL(0.C,A1,VIN2) 

TM    =    ICLIP    *    KT 

*  TNI    =    TM    -    TL 

TN2    =    TNI    *    (  1.0/EJ 

WM    =    PEALPL(0.0,A2,TN2) 

WMRP*    =    hf    *     (30. /3. 14159) 

kMRPPR   =     WMRP*/N 

VEMF    =    HP    *     K£ 

THETA    =    INTGRUO.CWM) 

THDEG    =    THETA    *     (180.0/3.14159265) 
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THCON    =    THRS7 
PwR    =    WM    *    TM 

************************************************** 

*  PROCEDURE    ICL1P    SIMULATES    PRESCENCE    GF  * 

*  FREEWHEELING    LICCE    BY    CLIPPING   THE    NEGATIVE  * 

*  GOING    CURRENT    AT    C    AMPS  * 
************************************************** 

PROCEDLRE    ICLIP    =    CLPCRUMJ 

IF    ilH    .LE.    G.Ci     ICLIP    =    0,0 

IF    UM    .GT.    0.0)     ICLIP    =    IW 
ENDPROCEDURE 


***************************************************** 

*  PROCEDURE    VCLF     ESTABLISHES    LIMIT    CYCLE    BEFAVICR.       * 

*  IF    THE    VELOCITY    IS    BELOW    THE    ESTABLISHED    SETPCINT,* 

*  THE     INPUT    VOLTAGE    IS    PULSED    ON    ( VCLP    =    30)    ANC    IF    * 

*  THE     VELOCITY     IS    AEJVE    VTCL    +    VCMC,    THE    INPUT  * 

*  VOLTAGE    IS    PLLSEO    OFF    (VCLP    =   O.OJ  * 
***************  ************************************** 


PROCEDLRE    VCLP    =    VVVUERR,     VTCLJ 
MVTCL   =    -l.C    *    VTOL 
IF     (VERR    .GT.    VTCL)     VCLP    =    20.0 
IF     (VERR    .LT.    MVTOLJ    VCLP    =    0.0 

ENDPROCEDURE 

************************************************** 

*  TNI    SETS    THE    SIGN   OF    THE    LOAD    TORQUE    TO    ENSURE   * 

*  THAT    THE    TQRCLE    OPPOSES    THE    MOTOR    TORCUE    WHEN      * 

*  THE    fCTOF    QPEMTES    IN    THE    FORWARD    DIRECTION  * 

*  AND    ADDS    TO     TFE    MCTOR    TORCUE      WHEN  * 

*  THE    /"OTOR    OPERATES    IN    THE    REVERSE   DIRECTION  * 
************************************************** 

PROCEDLRE    TM  =  f  H  EwC(  V  IN  ,T*  ,TL  J 
IF(VIN.LT.O.C)    GO   TO    10 
TNI    =    If    -    TL 
GO    TO    15 
10    TNI    =    TM    *    TL 
15    CGNT1NLE 
ENDPROCEDURE 

************************************************** 

*  PWRA    HAS     EEEN    ADOED    TO    THE    PROGRAM    TO    ACCCUNT      * 

*  FOR    THE       KONLUEAP    ASPECT    OF    THE    MODELLEC  * 

*  MOTOR'S    OLTPLT    POWER    FUNCTICIS  * 
************************************************** 


PROCEDLRE    PWRA=T«IDDL( PWR, TL.TLPJ 
IF(TL.GE.TLP)    GC    TO    20 
PWFA    =     PWR 
GO    TO    25 
20    TwID   =    4.5G    *     (TL    -    TL P J 
PWREX    =    EXP(-TWID) 
PWFA   =    PWR    *    PWREX 
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25    CONTINUE 
ENDPROCEDURE 

************************************************* 

*  THIS    PROCEDURE    *.AS    ADOEO    TO    RESET    THE    ROTGR'S    * 

*  POSITION    AFTER    IT   REACHES    360    DEGREES  * 

*  BACK    TO    0    DECREES  * 
****** *********** ******* *********************** ** 

PROCEDURE    THRST  =  RESET{  JFACTHDEGJ 
TS    =    JFAC   *    3*0.0 
THRS7    =    THDEG-    TS 
IF<ThRST.UT.36C.0J    GO    TO    40 
JF4C   =    JFAC    *    1  .0 
40    CONTINUE 
ENDPRGCEOURE 


TERMINAU 

TITUE    BASIC    CC    MOTOR    SYSTEM 
TIMER    FINTIM    =    .015,     OUTCEU    =    .00004, 
PRDEL    =    .000C4,    DELT    =    .000001 
METHCD    1RAPZ 

*  OUTPUT    WNRPM,    THETA,     THRST 
PRINT    ICUIP,^CUP,    dMRPM,     THRST 
UABEU    MOTCR    SFEEC    Do  E    TO    STEF    INPUT 
PAGE    MERGE 

*  PAGE    XYFUOT 

END 

STOP 
ENDJOB 
/* 
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APPENDIX  C 
DATA  ANALYSIS  PROGRAM 


$ENTR* 

$JOB  CLRR,XREF 

C 

C    THIS    PROGRAM    15    CESIGNED    TO    YIELD   THE    AVERAGE    CURRENT 

C    AND    THE    RMS    CLRRENT       AS    WELL    AS    K-FACTGR    FCR    GIVEN 

C    SIMULATICN    OF    A    COMPUTER    KOOEL    OF    A    PwM 

C    CONTROLLED    6FISHLESS    DC    MCTGF 

C 

FEAL    JO, ICC, RES,  INCU C  ,RPMAVE , PI , T 1ME , T IM E2, FR Eg 

REAL     EO,  IAV,I2RMS,A,B,ZN,KFAC,IRMS,ICN,ICFF,I/SIN 

FEAL    VIN,TC,P,TL,DC,PW,ICEN 

FI    =    3.1415S 
C 

C    RES     IS   ThE    RESISTANCE    OF    THE    COILS 

C    INDUC    IS    THE     INDUCTANCE    EETWEEN    THE     TWO    TERMINALS 
C    VIN     IS    THE    INFLT    VOLTAGE    CF    THE    MOTuR 
C 

RES    =   3.54C 

INDJC   =     .C013 

VIN    =   30. C 
C 

C      FREC    =    FREQUENCY.    OF    PwM    *«AV  E    SHAPES 
C      DC    REFERS    TC    CbTV   CYCLE    CR    • CN    TIME'     OF    THE     WAVE 
C 

C      TL    IS    THE    LCAC    TCRQUE    FCR    THE    GIVEN    RUN 
C       10     IS    THE    PEAK    CLRRENT    ACHIEVED 
C       100    IS    THE    MINIMUM    CURRENT    SEEN    BY    THE    MGTOR 
C       TO    IS    THE    AMCINT   OF     TIME    THAT    THE    CURRENT    WAVESHAPE 
C      IS    LESS    THAN    ZERC 
C 

FREQ    =    1H3. 

LX      =    .6*4 

11   =    9.6 

RPMAVE    -     1400 

IC   =    4.47 

ICC    =    3.42 

TO   =     .OOC 
C 

C       ZH    IS    THE    SYSTEM    ELECTRICAL    TIME   CONSTANT 
C 

2N    =    INCLC/PES 

FW   =    (1.0/FREC)    *    DC/100. 

TIME    =    PW 

TIME2   =    l.C/FREQ    -    PW 
C 

C    EO    IS    THE    BACK    EMF 
C 

EC   =     1.631    *    RPMAVE    *.0558   *PI/3C.O 
C 

C    IN    IS    THE    THECRETICAL    MAXIMUM    OF    THE    CURRENT 
C    IM    IS    THE   THECRETICAL    MINIMUM    OF    THE    CURRENT 
C    A    IS    THE     PULSE    CN   TIME 
C    B    IS    THE    PULSE    CFF    TIME 
C 
C 
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c 

c 


c 
c 


IW    =    EQ/FES 

IN    =     (VIN    -    ECJ/RES 

4    =    ZN    *    4LCGMIN    -    IGOJ/UN    -    IOJ 
E    =    ZN    *    ALGGiUO    ♦    If)/(IOO    ♦    IMJJ 


F    =    A    +    E    ♦    TC 

IAV    =    (  U*IN)    -     (E*IMJi     /    P 

IOEN       =     ZMUC    -    IOOJ*(JN    ♦     IMJJ/P 

12RMS   =     U4*IM*IN    +    B*IM*IM)    -    I  DEN 

IRKS    =    SCF7U2RMS) 

KFAC    =    IFNS/IAV 


WRITE    (  6  tlCJ 
10         FGRM/TdJj'STATISTICS    FCR    PwM    CONTROLLED    CC    MCTGR:*) 

WRITE    (6,20)FPEQ 
20         FGRMATi  IX, 'FREQUENC/    =    •  ,F8.2  ,'  H  Z  '  ,  J 

WRIT  E    <6,22)TL 
22         FORMAT!  LJi'lJDAO    TORQLE    =    %f-8.2i 

WRITE    (6,50)     IN,  IM 
50         FORMAT    <1>,«IN:     •fF10.3,«     IM:    «,F10.3J 

WRITE    (  6,  1CCJ    IAV,IRMS 
100       FCFM/Tll*  ,'IAV    =     '^8.^,/,'     IRMS    =    «,F8.4) 

WRITE    (  6f  iCOJ    KFAC 
200       FCRM4T    UX,«K   FACTOR    =    «.F8. 5) 

STOP 

ENC 
SENTRY 
SENTRY 
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APPENDIX  D 
SAMPLE  OUTPUT  FOR  CSMP  SIMULATION 


TIME 

KLIP 

VCLP 

bPRPP 

.3 

.0 

30.000 

.0 

4.0C300D-05 

.3*860 

30.000 

1.1220 

8.0COCOC-05 

.7244  6 

30.000 

4.4361 

1.2COOOC-04 

1.C677 

30.000 

9.8656 

1.6CC00C-04 

1.3984 

30.000 

17.334 

2.0C000C-0  4 

1.7168 

30.000 

26. 765 

2.4COOOD-04 

2.023  0 

30.000 

38.086 

2.8C00OC-O4 

2.3171 

30.000 

51.221 

3.2C000C-04 

2.599  3 

30.000 

66.099 

3.6COOOD-04 

2. £697 

30.000 

82.647 

4.0C300C-04 

3.1285 

30.000 

100.79 

4.4COOOC-04 

3.3759 

30.000 

120.47 

4.8C000C-04 

3.6119 

30.000 

141.61 

5.2C000C-04 

3.8369 

30.000 

164.13 

5.6CC0OD-O4 

4.C510 

30.000 

187.99 

6.0COOOC-04 

4.2544 

30.000 

213.10 

6.4COOOC-04 

4.447  1 

30.000 

239  .  40 

6.8C000D-04 

4.62  9  5 

30.000 

266.84 

7.2C000C-O4 

4.6018 

30.000 

295.25 

7.6C00OC-O4 

4.9641 

30.000 

224.  £6 

8.0C0COC-O4 

5.1165 

30.000 

255.22 

8.40000C-04 

5.2594 

30.000 

286.68 

8.8C0GOD-O4 

5.3929 

30.000 

418  .86 

9.2C0C0C-04 

5\5172 

30.000 

451. £2 

9.60000C-04 

5.632  5 

30.000 

485.49 

l.OCOGOD-03 

5.7391 

30.000 

519.84 

1.04000C-03 

5.  £370 

30.000 

554.79 

1.0£000C-03 

5.926  6 

30.000 

590.21 

1.12000C-03 

6.C080 

30.000 

626.24 

1.16000C-03 

6.C814 

30.000 

662.  £3 

1.2C00OD-Q3 

6.1471 

30.000 

699.74 

1.2400OC-O3 

6.2052 

30.000 

737. CI 

1.2£000C-03 

6.2560 

30.003 

774.61 

1.3200OC-03 

6.2996 

30.000 

812.48 

1.3  6300  C-0  3 

6.336  2 

30.000 

850.59 

1.4COGOD-03 

6.3662 

30.000 

£68.89 

1. 44000 C-03 

6.3895 

30.000 

927.25 

1  .4£000C-03 

6.4066 

30.000 

965.93 

1. 520000-03 

6.4175 

30.000 

100  4.6 

1.560C0C-03 

6.422  5 

30.000 

1043.3 

1.6COCOC-03 

6.4217 

30.300 

1062.0 

1.64COOC-03 

6.4154 

30.000 

1120.6 

1.6  8000  C-0  3 

6.403  7 

30.000 

1159.2 

1.720C0D-03 

6.3868 

30.000 

1197.7 

1.76000C-03 

6.3650 

30.000 

1236.1 

1.80000C-03 

6.3384 

3C.0C0 

1274.4 

1.8^0000-03 

6.3072 

30.000 

1312.4 

1.880C0C-Q3 

6.2716 

30.300 

1350.2 

1.92000C-03 

6.2317 

30.000 

1387.8 

1.9*000 C-03 

5.9513 

.0 

1424.7 

2.0C000C-03 

5.543  1 

.0 

1459.2 

2.0^0000-03 

5  .  143  8 

.0 

149  1.3 

2.08000C-03 

4.7534 

.0 

1520.9 

2. 12000  C-0  3 

4.3722 

.0 

1548.2 
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2.U0Q0C-03  4.C000  .0  1573.2 

2.2C000C-03  3. £369  .0  1596.0 

2. 240000-03  3.2829  .0  1616.5 

2.280C0C-03  2.9380  .0  1635.0 

2.32000C-03  2.6022  .0  1651.3 

2.360COC-03  2.2755  .0  l6t>5.7 

2.4G000C-03  1.9579  .0  1678.1 

2. 440000-03  1.6493  .0  1668.5 

2.48000G-03  1.3497  .0  1697.2 

2.52000C-03  1.0591  .0  1704.0 

2.560000-03  .77751  .0  1709.1 

2.6COOOC-03  .5C474  .0  1712.4 

2.64000C-0O  .24080  .0  1714.2 

2.660COC-03  .0  .0  1714.3 

2.72000C-03  .0  .0  1713.8 

2.7*0C0D-03  .0  .0  1713.3 

2.6C000C-03  .0  .0  1712.7 

2.84000C-03  .0  .0  1712.2 

2. 880000-03  .0  .0  1711.6 

2.92O0OC-O3  .0  .0  1711.1 

2.960COC-03  .0  .0  1710.5 

3.0COOOC-03  .0  .0  1709.9 

3.04000C-03  .0  .0  1709.3 

3.060COD-03  .0  .0  1706.7 

3.12000C-03  .0  .0  1708.1 

3.16000C-03  .0  .0  1707.5 

3.2COOOO-03  .0  .0  1706.9 

3.240C0C-03  ..0  .0  1706.3 

3.2SO0OC-O3  .0  .0  1705.6 

3.320C0C-03  .0  .0  1705.0 

3.36000C-03  .0  .0  1704.3 

3.4C0C0D-03  .0  .0  1703.7 

3.44000C-03  .0  .0  1703.0 

3.48000C-03  .0  .0  1702.3 

3.52C0OD-03  .0  .0  1701.6 

3.560C0C-03  .0  .0  1701.0 

3.6CC00O03  .0  .0  1700.3 

3.640C0C-03  .0  .0  1699.6 

3.68000003  .0  .0  1698.8 

3.72000D-03  .0  .0  1698.1 

3.760000-03  .0  .0  1697.4 

3.80000C-03  .0  .0  1696.7 

3.84C00D-03  .0  .0  1695.9 

3.8*0000-03  .0  .0  1695.2 

3.92000C-03  .0  .0  1694.4 

3.960COC-03  .0  .0  1693.6 

4.00000C-03  -0  .0  1692.9 

4.040C0D-03  .0  .0  1692.1 

4.08000C-03  .0  .0  1691.3 

4.12000C-03  .0  .0  1690.5 

4.16000D-03  .0  .0  1689.7 

4.2CO0OC-03  .0  .0  1688.9 

4.240C0C-03  .0  .0  1668.1 

4.260COD-03  .0  .0  1687.3 

4.32000C-03  .0  .0  1686.4 

4.360000-03  .0  .0  1685.6 

4.4C0C0C-03  .0  .0  1684.7 

4.44000C-03  .0  .0  1663.9 

4.4*0000-03  .0  .0  1683.0 

4.52000C-03  .0  .0  1682.1 

4.560COC-03  .0  .0  1681.3 

4.6C0CO0-O3  .0  .0  168C.4 

4.64000C-03  .0  .0  1679.5 

4.6*0000-03  .0  .0  lo78.6 

4.7iOOOC-03  .0  .0  1677.7 

4.76000C-03  .0  .0  1676.7 

4.8C0000-03  .0  .0  1675.8 
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4.84000C-03  .0  .0  1674.9 

4.860000-03  .0  .0  1673.9 

4.92000C-03  .0  .0  1673.0 

4.960C0D-03  .0  .0  1672.0 

5.0C000C-O3  .C  .0  1671.1 

5.04000C-03  .0  .0  167C.1 

5.080000-03  .0  .0  1669.1 

5.120COC-03  ..0  .0  1666.1 

5.1*CC0i>03  .0  .0  1667.1 

5.2CGC0C-03  .0  .0  1066.1 

5.24000C-03  .0  .0  1665.1 

5.26CC0D-03  .0  .0  1664.1 

5.32000D-03  .0  .0  1663.1 

5.36000C-03  .0  .0  1662.1 

5.4CCOOD-03  .0  .0  166  1.0 

5.44CC0C-03  ..0  .0  1660.0 

5.46CC0C-03  .0  .0  1658.9 

5.520C0C-03  .0  .0  1657.9 

5.56000C-03  .0  .0  1656.8 

5.6CCC0D-03  .0  .0  1655.7 

5.64CC0C-03  .0  .0  1654.6 

5.66000C-03  .0  .0  1653.5 

5.72C00D-03  .0  .0  1652.4 

5.760C0C-03  ..0  .0  1651.3 

5.8CC0OC-O3  .0  .0  1650.2 

5.84C00C-03  .0  .0  1649.1 

5. 860000-03  .0  .0  1647.9 

5.920000-03  .0  .0  1646.8 

5.96CC0C-03  .0  .0  1645.7 

6.0COCOC-03  .0  .0  1644.5 

6.C4000C-03  .0  .0  1643.3 

6.08000C-03  .0  .0  1642.2 

6.120C0D-03  .0  -0  1641.0 

6.U0COC-O3  .0  .0  1639.8 

6.20000C-03  .0  .0  1638.6 

O.24000D-03  .0  .0  1637.4 

6.26G00D-03  .,0  .0  1636.2 

6.32000C-03  .0  .0  1635.0 

6.360COC-03  .0  .0  1033.6 

6.4G000C-O3  .0  .0  1632.5 

6.44CC0D-03  .0  .0  1631.3 

6.46000C-03  .,0  .0  1630.1 

6.52000C-03  .0  .0  1628.8 

6.560000-03  .0  .0  1627.5 

6.6C0C0C-03  .0  .0  1626.3 

6.640C0C-03  .0  .0  1625.0 

6.66CC0C-03  .0  .0  lc23.7 

6.72000C-03  .0  .0  1622.4 

6.7*CC0D-03  .0  .0  1621.1 

6.8COOOC-03  .0  .0  161$. 8 

6.84000C-03  .0  .0  1618.5 

6.860000-03  .0  .0  1617.2 

6.92000C-03  .0  .0  1615.8 

6.9*0C0C-03  .0  .0  1614.5 

7.0COOOO-03  .0  .0  1613.2 

7.04000C-03  .0  .0  1611.8 

7.060000-03  .0  .0  lbl0.4 

7.12000C-03  .0  .0  1609.1 

7.16000C-03  .0  .0  1607.7 

7.2C00OD-O3  .0  .0  1606.3 

7.24000C-03  .0  .0  1604.9 

7.260COC-03  .0  .0  1603.5 

7.32000C-03  .0  .0  1602.1 

7.3c000C-03  .0  .0  160C.7 

7.4COOOO-03  .0  .0  1599.3 

7.44000C-03  .0  .0  1597.9 

7.46000C-03  .0  .0  1596.4 
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7.520000-03 

.0 

7. 560000-03 

.0 

7.6C0OOD-O3 

.0 

7.64CC0C-03 

.0 

7.66000C-03 

.0 

7.72COOC-03 

.0 

7.763000-03 

.0 

7.8COOOD-03 

.0 

7. 840000-03 

.0 

7. 860000-03 

.0 

7.92C00C-03 

.0 

7.96000C-03 

.0 

8.0C0COD-O3 

.0 

8.040000-03 

.0 

8.08000C-03 

.0 

8.120000-03 

.0 

8.163C0C-03 

•  0 

8.2C00OC-O3 

.0 

8.240000-03 

.0 

8.28000C-03 

.Q 

8.320000-03 

.0 

8. 360000-03 

.0 

8. 400000-03 

.0 

8.440000-03 

.0 

8. 460000-03 

.0 

8.52000C-03 

.0 

8. 560000-03 

.0 

8.60000C-03 

.0 

8.64CC0D-03 

.0 

8. 660000-03 

..0 

8.72000C-03 

.0 

8.760000-03 

.0 

8.8COOOC-03 

.0 

8.640000-03 

.0 

8.860000-03 

.0 

8.92000C-03 

..0 

8.96CC0D-03 

.0 

9.0COOOC-03 

.0 

9.04000C-03 

.0 

9.060000-03 

.0 

9.12000C-03 

.0 

9.160000-03 

.0 

9.2COCOC-03 

.0 

9.24000C-03 

.0 

9. 260000-03 

.0 

9. 320000-03 

.0 

9. 3 6000 0-0 3 

.0 

9.4COOOD-03 

.0 

9.44000C-03 

..c 

9.460000-03 

.0 

9.520COC-03 

.0 

9.56000C-03 

..0 

9.6COCOD-03 

.0 

9.640000-03 

.0 

9.68000C-03 

..0 

9.72000D-03 

.0 

9. 760000-03 

.0 

9.8C00OC-O3 

.0 

9. 840000-03 

.0 

9.86000C-03 

.0 

9. 920000-03 

.0 

9. 960000-03 

.0 

1.0C000C-02 

.Q 

1.0C400D-02 

..0 

1.0C8C0C-02 

.0 

1. 012000-02 

.0 

1.0  16000-02 

-.0 

.0 
.0 
.0 
.0 

•  0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 

:8 

.0 

•  0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 


.1 

.o 

.1 


1595.0 

1593.5 

1592 

159C 

15a9 

1587 

1586.1 

1584.6 

1583.1 
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